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We theoretically analyze the carrier capture and distribution among the available energy levels of a 
symmetric semiconductor quantum dot under continuous-wave excitation resonant with the barrier 
energy levels. At low temperature all the dot level-occupations but one decrease monotonically 
with energy. The uncovered exception, corresponding to the second (dark) energy level, displays a 
steady-state carrier density exceeding that of the lowest level more than a factor two. The root cause 
is not radiative recombination before relaxation to lower energy levels, but at the opposite, carrier 
trapping due to the symmetry-induced suppression of radiative recombination. Such a behaviour 
can be observed by collection-mode near-field optical microscopy. 



Understanding how electron-hole pairs distribute 
among the available energy levels of semiconductor 
nanostructures after optical excitation or current injec- 
tion is crucial for the development of novel quantum de- 
vices One tool frequently used to explore the energy 
distribution of electron-hole pairs is photoluminescence 
spectroscopy (PL). It provides a direct measurement of 
the optical density of states times the excitonic popula- 
tion density as a function of energy. Detailed simulations 
of Zimmermann et al 0, 0] have clarified many aspects 
of the intriguing nonequilibrium dynamics determining 
the distribution of excitonic populations among the avail- 
able energy levels in quantum wells (QWs) with interface 
roughness. Already at temperatures around 40-50 K the 
individual occupation of states is close to the Boltzmann 
distribution since all states can frequently emit and ab- 
sorb phonons to reach equilibrium before radiative emis- 
sion. At lower temperature more localized states display- 
ing smaller phonon scattering rates cannot equilibrate 
and their occupations deviate from the Boltzmann dis- 
tribution, displaying a reduced energy dependence. Far- 
field PL spectroscopy provide average measurements over 
a distribution of different emission sites. The opportu- 
nity for important insight is often lost by the inability 
to resolve finer details within this distribution. Scan- 
ning near-field optical microscopy (SNOM) combines the 
advantages of nanometric resolution of scanning-probe 
microscopy with the unique possibility of characteriz- 
ing quantum systems offered by optical spectroscopy. 
The SNOM ability to resolve the individual quantum 
constituents of semiconductor nanostructures has been 
widely demonstrated [1, H, If the spatial near-field 
resolution falls below the extension of confined quan- 
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tum systems, spatially resolved PL maps out the spa- 
tial probability distribution of the wave function times 
the corresponding level occupation 0, Q. The matrix 
elements governing the light-matter coupling are a con- 
volution of the quantum states with the exciting elec- 
tromagnetic field. This convolution implies that exciting 
a direct gap bulk semiconductor with a light field of a 
given wavevector resonant with the energy gap results 
in exciting cxcitons (i.e. bound state interband optical 
transitions) within the same wavevector. Succeeding in 
confining the optical excitation to a very small volume 
below the diffraction limit, implies the presence of opti- 
cal fields with high lateral spatial frequencies, determin- 
ing coupling matrix elements that can differ significantly 
from far-field ones [1, d, [H [HI, [H, [H . The most strik- 
ing manifestation of these effects is the breaking of the 
usual optical selection rules and the possibility to excite 
dark states whose optical excitation is forbidden by sym- 
metry in the far field. Spatial maps of dark states in 
semiconductor nanostructures were simulated for high- 
resolution SNOM in absorption and two-photon experi- 
ments 1^, [13, m, [3 • Moreover, dark states are not able 
to emit light in the far-field, for they generate only fields 
with high in-plane wavevectors corresponding to evanes- 
cent waves. Hence their lifetimes at low temperatures 
result significantly increased with respect to their bright 
counterparts. 

In this letter we theoretically investigate the impact of 
the reduced dark-states relaxation on the distribution of 
electron-hole pairs among the available energy levels after 
a far-field continuous wave excitation. We find that this 
reduced relaxation significantly throws off balance the PL 
population dynamics resulting in striking deviations from 
equilibrium. In particular we find that, even at steady 
state, the occupation of the dark first excited state of a 
symmetric artificial atom can exceed that of the ground 
state. Usually larger populations for higher energy levels 
in nanosystems are observed as consequence of bottleneck 
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FIG. 1: (color online), (a) Near-field PL signal (collection- 
mode) as a function of photon energy and beam position (line 
scan shown in the inset) obtained under uniform illumination 
of the sample at the energy of the QW Is exciton. T=2 K, 
FWHM =40 nm. (b) Total absorption under local illumina- 
tion, along the same line, spatial resolution and temperature 
as in (a). 

effects due to efficient radiative recombination occurring 
before relaxation to the lowest energy-levels. Here the 
inverted occupations origin from the opposite effect: car- 
rier trapping due to the presence of nonradiative states. 
While the observation of such effects by means of far- 
field spectroscopy is prevented (dark states are not able 
to emit), a SNOM tip can collect the evanescent waves 
generated by the occupation of dark states. Hence near- 
field PL reveals to be an excellent laboratory to address 
general questions regarding nanoscale energy transfer in 
open quantum systems. 

In this work we consider quantum states confined 
in a naturally occurring quantum dot, or terrace, in- 
duced by monolayer fluctuations in the thickness of a 
GaAs/AlGaAs semiconductor quantum well. We adopt 
the usual envelope function and assume that the electron- 
hole wave functions factorizc into a center of mass and a 
relative part given by the Is state of quantum well cxci- 
tons Q. We assume that the radiative decay of excitons 
is not drastically altered by the presence of the SNOM 
tip. As pointed out in Ref. [3, this is a reasonable as- 
sumption since the photons can be emitted into any solid- 
angle direction, and the slightly modified photon density 
of states in the presence of the SNOM tip is not expected 
to be ofgreat importance. Within these assumptions, re- 
cently Q it has been shown that the near-field spectrally 
resolved PL signal collected by the tip can be written as 

IPL 

a 

(1) 

where o°"* = J (PR Eout{'R-)4'a (R) is the overlap of the 
exciton wavefunctions with the signal mode Eo^t (R-) sup- 
ported by the tip centered at the position Rout (collec- 
tion mode), nCa{uj) = T/[{uj~LJa)^ +r^]i where F is the 
dephasing rate of the exciton due to radiative emission 
and phonon scattering. In Eq. (1) Nq, are the diagonal 
terms of the exciton density matrix. The kinetic equation 
for the diagonal terms of the exciton density matrix can 
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FIG. 2: (color online), (a) PL spectra obtained in collec- 
tion mode {FWHM = 40 nm and T = 4 K). The spectra 
have been obtained centering the tip at one of the two emis- 
sion maxima of the first excited state (LS) and at the emission 
maximum of the ground state (C); (b) Total absorption under 
local illumination obtained with spatial resolution, at temper- 
ature and centering the snom tip as in Fig. 2(a); (c) PL spec- 
tra as in Fig. 2(a) where the radiative decay rates of all the 
levels are set two orders of magnitude lower than the lowest 
nonradiative decay rate, (d) Level occupations (normalized 
at the maximum) calculated for two different temperatures. 

be derived from the Hcisenbcrg equations of motion for 
the exciton operators under the influence of the system 
Hamiltonian Q . The relevant Hamiltonian is constituted 
by the bare electronic Hamiltonian of the semiconductor 
system, the interaction Hamiltonian of the semiconduc- 
tor with the light field and the interaction Hamiltonian of 
excitons with the phonon bath. In the following we con- 
sider a low-excitation regime, and according to the dy- 
namics controlled truncation scheme (l7| we include only 
states with one electron- hole pair (excitons). In addition, 
we neglect the possible coherent phonon states and the 
memory effects induced by the photon and phonon fields. 
The resulting set of kinetic equations is Q 

where 7/3^a are the phonon- assisted scattering rates. 
2Ta ~ ra+ X]/3 ip^a is thc total out-scattering rate, be- 
ing the rate for spontaneous emission, proportional to 
thc exciton oscillator strength: = ?'o |/ (PRipa (R)| ■ 
In this equation thc generation term depends on the spa- 
tial overlap between the illuminating beam and the exci- 

ton wavefunctions: Gq = ?'o |o™| Cai^^in) and o™, anal- 
ogously to 0°"* , contains the overlap of the exciton wave- 
functions with the signal mode £'i„(R) delivered by the 
tip centered at Ri„. In collection mode (far- field illumi- 
nation) the input field can be considered as uniform over 
thc QW plane: i?™(R) ^ E°„ . 

We apply the above developed theoretical scheme to 
calculate the individual occupations of exciton states con- 
fined in the dot after a continuous-wave far-field optical 
excitation resonant with the energy of the Is QW exci- 
ton (the dot barrier). The obtained occupations are then 
used to study the spatially and spectrally resolved (col- 
lection mode) light emission from thc dot. The effective 
potential felt by Is excitons used in our simulations is 
represented by a sample of (240 x 240) nm with a pro- 
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FIG. 3: (color online), ((color online). PL spectral line-scans 
calculated at three different temperatures, with FWHM =40 
nm. 

totypical interface- fluctuation confinement of rectangular 
shape with dimensions (60 x 90) nm, and monolayer fluc- 
tuations giving rise to a 6 meV effective confinement po- 
tential. The upper panel of Fig. 1 displays the near-field 
PL signal as a function of photon energy and beam posi- 
tion obtained after uniform illumination of the sample at 
the energy of the Is QW exciton in the absence of inter- 
face fluctuations, fixed as zero of energy {ivj ~ meV). 
The specific line scan is indicated in the inset. A Gaus- 
sian profile with FWHM = 40 nm of the electromagnetic- 
mode supported by the collecting tip has been assumed. 
The line scan clearly evidences the first excited state of 
the dot which is dark under far-field collection. We ob- 
serve that its spectral line is more intense and narrow 
(owing to the absence of radiative decay) than that of 
the ground state. The calculated PL spectra shows that 
dark-states can be observed by high-resolution SNOM 
in the usual collection mode configuration after nonres- 
onant far-field excitation, without the need of nonlinear 
optical interactions. Analogous results are expected for 
locally collected electroluminescence. For comparison we 
plotted in the lower panel the total absorption under lo- 
cal illumination (FWHM = 40 nm). As near-field PL 
is proportional to this quantity times the level occupa- 
tions, the comparison provides interesting information 
about the steady-state exciton populations. The two 
panels clearly indicate that the level occupation of the 
dark-state is significantly larger than that of the ground 
state, at the opposite to what predicted by the Boltz- 
mann distribution. The observed inverted occupations 
origin from symmetry-suppression of radiative decay of 
the dark-state and the quite small nonradiative scatter- 
ing at low temperature {T — 2K) for states well con- 
fined in the dot. This behaviour is better evidenced in 
panels 2(a) and 2(b). The spectra have been obtained 
centering the tip at one of the two emission maxima of 
the first excited state (LS) and at the emission maxi- 
mum of the ground state (C). In order to better spec- 
ify the origin of the inverted occupations, we calculated 
the PL spectra as in Fig. 2(a) except that we artificially 



FIG. 4: (color online). PL spectral line-scans calculated at 
three different spatial resolutions, at the temperature of 4 K. 

set the radiative decay rates of all the levels two orders 
of magnitude lower than the lowest nonradiative decay 
rate. Panel 2(c) displays the results showing that, as 
expected in this case, equilibrium is recovered and the 
near-field emission from the dot lowest energy level dom- 
inates even when the tip is centered on the maximum 
of the dark-state emission. Fig. 2(d) shows the calcu- 
lated level occupations at T = 2K and T = 30K. At 
T = 2K all the occupations with the striking exception of 
the second energy level decrease monotonically with en- 
ergy. This level displaying an occupation density which 
is more then a factor two that of the lowest energy level. 
At T = 30A' the monotonous behaviour is recovered for 
all the energy levels. We notice that the sum of the occu- 
pations obtained at T = 30Ar is significantly larger then 
that at T = 2A'. At low temperature a large fraction 
of the resonantly generated QW excitons (with a quite 
large radiative decay rate) recombine by radiative emis- 
sion, while at higher temperature phonon scattering low- 
ers this effect, increasing the carrier density captured by 
the dot. Fig. 3 shows the near-field {FWHM ^ 40 nm) 
PL spectral line-scans obtained at three different tem- 
peratures. Increasing the temperature all states are able 
to better thermalize, frequently emitting and absorbing 
phonons before radiative emission. The influence of spa- 
tial resolution on the PL spectral line-scans is shown in 
Fig. 4. This figure displays the line-scans calculated at 
T = 2K for three different spatial resolutions (indicated 
in the figure). Lowering spatial resolution, the signal 
from the dark state tend to disappear owing to cancel- 
lation effects in the matrix element o°"* governing the 
light-matter coupling. Finally, we test the influence of 
random interface fluctuations with a correlation length 
of the order of the exciton Bohr-radius 0, 0] . The total 
effective potential felt by excitons is obtained adding to 
the dot potential used for all the previous calculations a 
contribution modeled as a zero mean. Gauss distributed 
and spatially correlated process with a correlation length 
a = 8 nm and a width of the energy distribution vq ~ 0.6 
meV (about the 10% of the dot potential- well). Fig. 5(a) 
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FIG. 5: (color online), (a) Specific realization of the effective 
disordered potential used for the calculation of the PL images 
and spectra show, (b) PL energy-integrated image obtained 
after uniform illumination of the sample at the barrier energy 
and collecting locally {FWHM = 40 nm, T = 4 K). (c) PL 
spectra calculated centering the tip in the center of the dot 
and in the position indicated by a circle in the first panel 
{FWHM = 40 nm, r = 4 K). (d)Far-field PL spectrum 
{FWHM = 40 nm, r = 4 K). 

displays a specific realization of this potential. Fig. 5(b) 
displays the energy-integrated PL image {FWHM = 40 



nm, T = 2K) for the potential realization shown in panel 
5(a). It evidences the efficiency of the dot capture and 
near-field emission. Fig. 5(c) shows that the presence 
of the symmetry-breaking disordered potential slightly 
increases the linewidth of the second energy level, deter- 
mining a slight equilibration of the emission lines. Fig. 
5(d) displays the far-field PL spectrum. A very small 
peak (not present in the absence of disorder) appears at 
the energy of the second level, which thus remains mainly 
dark, demonstrating the robustness of dark-states with 
respect to this kind of disorder. Results obtained for 
other specific realizations of the random potential (here 
not shown), don't display qualitative differences. 

In conclusion, we have investigated the impact of the 
reduced dark-states relaxation on the distribution of 
electron-hole pairs among the available energy levels af- 
ter a far-field continuous wave excitation. The calculated 
near-field luminescence properties of these states depend 
critically on tip position, temperature and spatial reso- 
lution and clearly indicate the potentiality of near-field 
PL for adressing general questions regarding nanoscale 
energy transfer in open nanosystems. 
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